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амотлорское месторождение по своим запасам отно- 

сится к уникальным. Промышленная нефтегазонос- 

ность месторождения связана с терригенными отло- 
жениями верхней юры (пласты ЮВ,!, ЮВ, васюганской 
свиты), нижнего мела (пласты БВ. 22, БВ,» БВ БВ17-1з, БВ 
БВ» Вз, БВз) мегионской свиты; пласты БВ., БВ., БВЗ, БВ», 
БВ;, БВ, АВ,, АВ-, АВь, АВ, 5, АВ,  ванденской свиты; АВ 3, 
АВ," алымской свиты; пласт ПК; нижней части покурской 
свиты) и верхнего мела (пласты ПК, ;!- ПК, средней и верхней 
части покурской свиты). 

Продуктивная толща месторождения характеризуется высо- 
кой литологической неоднородностью. В пределах месторож- 
дения развиты различные типы коллекторов с контрастными 
свойствами. Пласт АВ! («рябчик») представлен слоистыми 
аргиллитами и глинистыми алевролитами, биотурбированны- 
ми песчаными алевролитами и глинистыми песчаниками, от- 
личается линзовидностью строения, высокой анизотропией 
свойств и расчлененностью, низкой проницаемостью ~ 
(8-10)-103 мкм?. Пласт БВ, образован шельфовыми отложе- 
ниями с высокой песчанистостью и средней проницаемостью 
0,5 мкм?. Пласты БВ; и БВ, 6.22 имеют клиноформенное строе- 
ние и представляют собой изолированные песчаные ловушки с 
разными уровнями водонефтяного контакта (ВНК). Пласт ЮВ, 
сложен переслаивающимися аргиллитами, песчаниками и 
алевролитами со средней проницаемостью 0,02 мкм?, Пласты 
АВ; 5 и ПК имеют газовые шапки. 

_ Самотлорское месторождение разрабатывается с 1969 г, за 


стовых систем. По состоянию на 01.01.10 г объектами разра- 
ботки являются пласты АВ,!, АВ,?, АВ, АВ, ;, АВ в, БВо- 
БВ,!, БВр, БВ, БВ1б.22 и ЮВ.. Пласты группы ПК в настоящее 
время не разрабатываются. 

Проектный коэффициент извлечения нефти (КИН) составляет 
в среднем 0,497: от 0,325 для пласта АВ, 1"? до 0,658 для пласта 
БВ,. Текущий КИН по месторождению в целом равен 0,352 при 
обводненности продукции 93,1 %. По основным объектам разра- 
ботки текущий КИН варьирует в диапазоне от 0,07 (объект АВ\!?) 
до 0,569 (объект БВр). На месторождении планируется провести 
91 600 скважино-операций, в том числе 18 458 гидроразрывов 
пласта (ГРП), 2 015 зарезок боковых стволов (БС). 

Направлениями совершенствования реализуемых систем раз- 
работки по эксплуатационным объектам являются формирова- 
ние элементов блочно-замкнутых систем, дополнительное ли- 
нейное разрезание с образованием однорядных систем, созда- 
ние элементов очагового и площадного заводнения в зонах пре- 
рывистых и сильно прерывистых коллекторов, оптимизация 
плотности сетки скважин, интенсификация отборов, осуществ- 
ление мероприятий по регулированию разработки, увеличение 
степени использования пробуренного фонда скважин. 

Для повышения эффективности разработки созданы и актив- 
но используются геолого-гидродинамические модели всех объ- 
ектов месторождения. Параметры расчетных сеток основных 
моделей приведены в табл. 1. 

Очевидно, что создание моделей такого месторождения яв- 
ляется сложной задачей. Выделим некоторые проблемы модели- 
рования, обусловленные особенностями Самотлорского место- 
рождения. 

1. Большие размеры залежи (50х80 км), что приводит к уни- 
кально большим размерам моделей. 
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Таблица 1 
Размерность геологической модели 


ячеек 


661 | 1071 240696540 


Размерность 
гидродинамической модели 


была построена трехмерная трехфазная гидро- 

динамическая модель с использованием про- 

граммного комплекса (ПК) Есйрве. 
Традиционным подходом к уменьшению 


числа расчетных ячеек модели является ремас- 
7768200 


638 | 1054 257 176182424 


2689808 штабирование исходной геологической модели. 


Е Е 662 | 734 | 297 | 14431467 


858880 Важно соблюдать баланс между практической 


75 | 75 | 531 | 530 | 209 | 58818870 


909480 | пользой и физической сущностью модели. Чрез- 


75 | 75 | 675 | 1062 235 168459750 


2849280 


50 | 50 | 396 | 365 | 107 | 15465780 


2. Высокая плотность сетки скважин, большое число горизон- 
тальных БС. Так, для создания геологической модели самого 
крупного объекта АВ,.; использована информация по 18 тыс. 
скважин, а в гидродинамическую модель загружены данные по 
истории добычи почти 14 тыс. добывающих скважин. В более 
чем половине фонда скважин осуществляется совместная экс- 
плуатация двух и более объектов. Такой объем промысловой ин- 
формации в модели с учетом 40-летней истории разработки мес- 
торождения не позволяет в настоящее время использовать про- 
граммные средства автоматической адаптации модели даже при 
наличии значительных вычислительных мощностей. 

5. Сложность создания модели насыщенности, так как пласты 
имеют гетерогенное строение. Например, залежь пласта АВ! 
является недонасыщенной нефтью, толщина переходной зоны 
в краевых областях достигает 40-50 м, при этом модель должна 
учитывать различие фильтрационных характеристик, высокую 
начальную обводненность продукции скважин и быстрые про- 
рывы воды без выработки геологических запасов на участках 
краевых зон. 

4. Большая неопределенность в оценке текущих пластовых 
давлений и уровня газонефтяного контакта (ТНК), свойств 
флюидов, начальных нефтенасыщенностей, насыщенностей 
коллекторов газовой шапки нефтью, отборов газа скважинами, 
свойств глинистых перемычек и других характеристик пластов и 
скважин. 

5. Сложность отражения в модели основных процессов, свя- 
занных с наличием газовых шапок: продвижение нефти в газо- 
вую шапку и возможное поднятие уровня ГНК в процессе раз- 
работки, так как газовая шапка пластов АВ\.5 занимает около 
130 км?, или примерно 25 % всей площади залежи, и влияет на 
выработку запасов. 

6. Сложность воспроизведения в модели различных геолого- 
технологических мероприятий (ТТМ), особенно ГРП. При вы- 
полнении ГРП в пласте АВ, 1? часто нарушается глинистая пере- 
мычка и подключается нижележащий пласт АВ? с более высоки- 
ми коллекторскими свойствами. 

7. Трудности в оценке объемов непроизводительной закачки 
воды нагнетательными скважинами при моделировании, так как 
при эксплуатации нагнетательных скважин возможно образова- 
ние техногенных трещин, в скважинах распространены заколон- 
ные перетоки. 

Рассмотрим подходы к решению указанных проблем на при- 
мере геолого-гидродинамической модели пластов АВ 5, ОДНОГО 
из основных объектов Самотлорского месторождения. 3р гео- 
логическая модель была создана в 2007 г для пересчета геологи- 
ческих запасов нефти и растворенного газа. На ее основе в рам- 
ках работы «ТЭО КИН Самотлорского месторождения» в 2008 г 


мерное упрощение (загрубление модели) может 
привести к полной потере неоднородности и не 
позволит получить адекватное отражение про- 
цессов фильтрации в пласте. Для сохранения исходной неодно- 
родности вертикальное укрупнение модели выполнялось нерав- 
номерно, по специальному алгоритму, основанному на исполь- 
зовании геолого-статистических разрезов (ГСР) и индивидуаль- 
ных характеристик слоев исходной геологической сетки. В ре- 
зультате средняя толщина ячеек в целом в модели составила 3,5 м, 
а в нефтенасыщенной части ~ 2,5 м. Это позволило уменьшить 
число активных ячеек в модели до 4,7 млн. и осуществлять «про- 
гон» истории за 10-12 ч на 20 процессорах Ваае. 

Огромное число скважин затрудняет процесс настройки. Так, 
для адаптации модели по забойным давлениям последнего исто- 
рического периода (2003-2009 гг) было обработано 378 тыс. 
определений забойных давлений, из них использовано в работе 
80935. Поэтому адаптация гидродинамической модели выполня- 
лась по группам скважин. За основу формирования групп была 
взята схема кустования скважин, выделено 260 групп по 50-60 
скважин в каждой. При настройке модели использовалась техно- 
логия укрупнения ячеек Соагѕеп. Модель была разделена на три 
примерно равнозначных блока с поочередной настройкой каж- 
дого из них. При этом на двух других блоках размеры ячеек по 
латерали многократно укрупнялись (рис. 1). 


936936 


Рис. 1. Схематическое расположение групп скважин и блоков, по кото- 
рым выполнялась адаптация геолого-гидродинамической модели объ- 
екта АВ. _5 
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Такой подход позволяет избежать неопределенностей, возни- 
кающих при фиксации перетоков между ячейками в неадаптиро- 
ванной модели (как это происходит при разрезании модели 
крупного объекта на секторы и поочередной их настройки), и 
сократить расчетное время вследствие кратного уменьшения 
размеров модели. 

Для учета неоднородного характера насыщения на этапе гео- 
логического моделирования было выделено по четыре типа кол- 
лекторов с четырьмя типами петрофизических зависимостей для 
пласта АВ! и пластов АВ, 5 (рис. 2). 


Проблемы неопределенности данных имеют огромное значе- 
ние. Например, наличие газовой шапки в пластах АВ! -5 свиде- 
тельствует о близких величинах начального пластового давления 
и давления насыщения нефти газом. Однако по результатам ана- 
лизов глубинных проб нефти давление насыщения примерно на 
40 % ниже начального пластового давления на ГНК. Для вос- 
произведения начальных условий залежи были выполнены ре- 
комбинации устьевых проб нефти, что позволило построить со- 
гласованную РУТ-модель с распределением давления насыщения 
в зависимости от глубины. 
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Рис. 2. Петрофизические зависимости нефтегазонасыщенности от превышения над зеркалом воды для пластов АВ. 3, АВ, з, АВ, (а) и АВ. 1-2 (6) 


Самотлорского месторождения 


При переходе к фильтрационной модели кубы начальной 
нефте- и газонасыщенности приведены в равновесное состоя- 
ние, схема инициализации показана на рис. 3. 

При значительном укрупнении исходной геологической сетки 
необходимо вносить поправки в функции относительных фазо- 
вых проницаемостей (ОФП), так как влияние численной диспер- 
сии может привести к серьезным погрешностям в результатах 
расчетов. В рассматриваемом случае была определена погреш- 
ность расчета КИН, которая составила 3 %, в исходные кривые 
ОФП внесена соответствующая поправка. 

Для более детального учета особенностей коллекторов с раз- 
личными фильтрационными свойствами в модели выделено пять 
регионов насыщенности, использовалась опция масштабирова- 
ния критических значений водо- и нефтенасыщенности в зави- 
симости от проницаемости. 


Для оценки вертикальной анизотропии проницаемости про- 
веден ряд вариантных расчетов, показавших, что динамика и 
уровень вычисленных пластовых давлений воспроизводится в 
модели близко к фактическим при соотношении вертикальной и 
горизонтальной проницаемостей по верхнему пласту 1:500, по 
нижним пластам - 1:100 и 1:10. 

Одной из основных задач модели является адекватное отобра- 
жение процессов, происходящих в пластах. В связи с отсутстви- 
ем на месторождении достоверных данных по определению 
коэффициента остаточной нефтенасыщенности К; в коллекто- 
рах газовой шапки возможные потери нефти от внедрения в га- 
зовую шапку оценивались при помощи имитационных многова- 
риантных расчетов с различными значениями К, В газовой 
шапке. Часть задач, например оценка влияния разгазирования 
нефти на снижение дебита, была решена при помощи имита- 


ционных моделей небольших размеров, позволяющих 
исследовать эффекты на детальной сетке (10-20 м). 
Имитация ГТМ в модели имеет условный характер. 
При выбранном размере ячеек (150х150 м) приходится 
делать значительно больше допущений, поэтому нельзя 
утверждать, что выбранные подходы можно применять 
для достоверных прогнозов по каждой отдельной сква- 
| жине. Моделирование ГРП и других ГТМ осуществля- 
лось изменением скин-фактора. Как уже отмечалось, 
особую сложность представляло воспроизведение эф- 
фекта, получаемого при осуществлении ГРИ в пласте 


зации модели пластов АВ. ; 


АВ, "2, когда разрушается глинистая перемычка и вклю- 
чается в работу нижележащий пласт АВ). Была пос- 
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троена карта глинистой перемычки между указанными пластами, 
и с использованием обратной зависимости от толщины этой пе- 
ремычки подключалась (или не подключалась) соответствующая 
часть нижележащего пласта. 

Непроизводительная закачка была оценена путем расчета по 
методу материального баланса, который показал, что при нуле- 
вой активности законтурной водоносной области объемы закач- 
ки, зафиксированные в отчетности, на 8 % превышают физичес- 
кие возможности залежей. Для упруговодонапорного режима с 
активной законтурной водоносной областью объем непроизво- 
дительной закачки по пластам АВ ; составляет 20-30 %, в связи с 
чем было решено не воспроизводить в модели фактическую за- 
качку в полном объеме, установив для нагнетательных скважин 
ограничения по давлению. Все это позволило выполнить адапта- 
цию основных параметров модели объекта АВ.; в целом и по 
группам скважин с погрешностью 5-10 % (рис. 4). 


определения текущего характера насыщения. Результаты его ана- 
лиза подтвердили полученную в модели динамику насыщенности 
(табл. 2). Расхождение текущей насыщенности пластов по резуль- 
татам моделирования и по данным анализа керна составило 3-15 %. 

Однако возможности использования моделей с ячейками 
150х150х3,5 м для решения многих задач, таких как оптимиза- 
ция заводнения, формирование адресных программ ГТМ, обос- 
нование бурения, особенно при уплотненной сетке скважин, и 
других – ограничены. Это обусловлено существенным загрубле- 
нием расчетных сеток при переходе от геологической модели к 
гидродинамической. С другой стороны, использование деталь- 
ной сетки в полномасштабных моделях приводит к лавинооб- 
разному росту числа ячеек, неустойчивости расчетов в мелких 
ячейках и необходимости использовать многопроцессорные и 
многоядерные вычислительные кластеры для выполнения дли- 
тельных вариантных расчетов. 
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Рис. 4. Результаты адаптации модели объекта АВ. .; к фактической истории по текущей добыче нефти (а) и добыче нефти по группам скважин (6) 


Учитывая большие неопределенности при распреде- Таблица 2 


лении в отчетности добычи между пластами, адапта- 
цию по пластам не проводили. В данном случае имен- 
но модель является инструментом, позволяющим более 
полноценно и разносторонне решить эту проблему. 


объект 


Исследуемый 


модель 


Текущее насыщение в среднем 

по модели, % (оценка 2009 г., исследованной скважине, % 
до отбора керна) (оценка 2010 г.) 

| нефтью | газом | водой | Всего | нефтью | газом | водой | 

3 35 62 100 - 38,9 61,1 


Текущее насыщение в 


Всего 


Полномасштабные геолого-гидродинамические керн 


24,3 73,8 


модели Самотлорского месторождения используются 


модель 


26 65 100 35 65,1 


керн 


30,4 | 68,2 


при решении актуальных задач прогноза добычи 


модель 


26 62 100 24,5 75,5 


нефти, оценки перспектив добычи газа из газовых керн 


шапок, влияния этих процессов на выработку запасов 


модель 


66 100 


13,7 


61,3 


нефти и др. Так, по результатам моделирования было Керн 


247 


59,9 


отмечено поднятие уровня ГНК на 10-30 м по пластам 


модель 


63 100 28 


65,7 


керн 
АВ, з, локализованы и определены текущие запасы 


нефти. В настоящее время пробурено несколько скважин, под- 
твердивших наличие нефти выше первоначального уровня ГНК. 

Результаты моделирования позволили объяснить газопроявле- 
ния в скважинах, перфорированных в чисто нефтяной зоне 
пластов АВ, з и имеющих глинистые перемычки, отделяющие 
объект АВ, , от вышележащих пластов АВ, и АВ, "?. Причиной 
может служить вертикальная миграция газа в нижележащие плас - 
ты через локальные «окна слияния» с последующим его проры- 
вом в интервал перфорации скважин. 

Достоверность модели подтверждена исследованиями керна. 
Так, в 2009-2010 гг. в одной из скважин был отобран керн для 


25,2 


67,4 


Основным преимуществом секторных моделей является то, 
ЧТО ОНИ могут быть сколь угодно детальными и использоваться 
для работы с каждой отдельно взятой скважиной. Для ряда участ- 
ков Самотлорского месторождения в последние годы были соз- 
даны секторные модели, однако, к сожалению, единая концепция 
создания секторных моделей для Самотлорского месторождения 
пока не сформирована. 

Однако использование секторных моделей имеет некоторые 
ограничения. 

1. Сложность создания секторных моделей для пластов с газо- 
выми шапками. 
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2. Нарушение целостности моделируемого объекта, в результа- 
те чего задание краевых условий на границах моделей может 
быть некорректным, что приводит к ошибкам моделирования. 

3. Накопленные в секторных моделях результаты исследова- 
ний скважин и межскважинных пространств не передаются в 
полные модели. 

4. Секторные модели создаются для решения узконаправлен- 
ных задач и имеют «короткую жизнь». 

В связи с этим построение полномасштабных моделей по- 
прежнему имеет большое значение даже в случае, если система 
секторных моделей плотно покрывает всю залежь. Для обоснова- 
ния и выбора вариантов разработки, составления проектных тех- 
нологических документов, формирования карт распределения 
текущих запасов, изобар и решения других задач необходимо 
создавать полномасштабные модели. 

Для Самотлорского месторождения актуально создание пол- 
номасштабной модели всех объектов разработки с детальностью 
секторных моделей. Число расчетных ячеек в такой модели 
может составлять до 107-109, что обусловлено крупными разме- 
рами залежей, большим числом объектов разработки и скважин, 
эксплуатирующих несколько объектов совместно, а также жела- 
нием создавать модели с минимальным загрублением исходных 
геологических моделей. 


Возможным решением проблемы полномасштабного и при 
этом детального моделирования крупных пластовых систем яв- 
ляется применение технологии раздельно-одновременного сек- 
торного моделирования". Она основана на декомпозиции исход- 
ной геологической модели на несколько секторных геолого-гид- 
родинамических моделей с последующим сопряжением сектор- 
ных моделей по потокам флюидов, полям давлений и насыщен- 
ностям, что обеспечивает целостность моделируемого объекта. 
Сопряжение секторных моделей проводится итерационно на каж- 
дом временном шаге или через определенный интервал по специ- 
альному алгоритму. В настоящее время эта технология апробиру- 
ется на модели пластов АВ, < Самотлорского месторождения. 

В ближайшие годы (2011-2014 тг) перед ООО «ТННЦ» стоит 
задача последовательного обновления, с учетом накопленного 
объема информации, петрофизической, геологической и гидро- 
динамической моделей. Успешная апробация новой технологии 
позволит создать модель с минимальным загрублением геологи- 
ческой сетки, снизить требования к ресурсам компьютеров, рас- 
параллелить работы в проектной группе и повысить качество 
адаптации. Наличие единой постоянно действующей полномас- 
штабной модели такой степени детальности, пригодной для ре- 
шения любых задач, обеспечит эффективное сопровождение 
разработки Самотлорского месторождения. 


1 Принцип итерационного сопряжения секторных моделей для полномасштабного моделирования больших и гигантских пластовых систем/А.С. Бордзиловский, 
И.С. Игнатов, С.В. Костюченко, Е.И. Шапиева//Нефтяное хозяйство. - 2009. – № 8. - С. 42-46. 
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